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Abstract 

In the previous work, we investigated the rotation curves of galaxies assuming 
that the dark matter consists of ultra light boson appearing in 'late time 
phase transition' theory. Generalizing this work, we consider the effect of 
visible matter and classify the types of rotation curves as we vary the fraction 
of the mass and extention of visible matter. We show that visible matter, 
in galaxies with flat rotation curves, has mass fraction 2% ~ 10% and it is 
confined within the distance fraction 10% ~ 20%. 
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I. INTRODUCTION 



Recently, motivated by the large scale structure, an ultra light Nambu-Goldstone boson 
was introduced in dark matter physics under the name of "late time cosmological phase 
transition" ffl]. The basic idea of the theory is that if a phase transition happened after 
decoupling, one can avoid the constraint imposed by the isotropy of the microwave back 
ground. 

Subsequently, one of the authors of this paper tried to apply this idea to the rotation 
curve of galactic haloes. If the phase transition occurred so late, the universe was already 
big, therefore the relevant particle must be such that (i) its Compton wave length provides 
the scale of interest since there is no other length scale of this kind, (ii) it should be a major 
component of dark matter || to be responsible for the structure formation. 

If the dark matter consists of this particle whose mass is, say, m ~ 10~ 24 eV and if the 
dark matter density in the galaxy is about 10~ 25 g/cm 3 , then the inter particle distance is 
of order 10~ 13 cm, while the Compton wave length is of lOpc order. Therefore the wave 
functions of the particles are entirely overlapping and a galactic halo is a giant system of 
a bose liquid. In this context, one must consider the dark matter distribution quantum 
mechanically and this is done in ref. 0. 

However, there was a point that had to be clarified further in that work. All the rotation 
curves obtained there were slightly increasing. While this is consistent with most of the 
data f|], H, there are a few galaxies whose rotation curves are decreasing. In fact, in [pf, 
the effect of the visible matter was not properly considered. Although most of the mass of 
a galaxy is attributed to the dark matter, the shape of the rotation curve can depend on 
the mass profile of the visible matter, especially near the core of a galaxy. In this paper, we 
consider the effect of the visible matter to the shape of the rotation curve more carefully in 
the scheme of the ref. 0. 

We briefly recapitulate the main idea of the paper ||. The dark matter that is consist 
of ultra-light boson is in a condensation state. Usually bosons condense in the ground state, 
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because in atomic physics situation, matter is always coupled to electromagnetism (cooling) 
and this efficiently lower the energy of a system. In our case, however, it was shown |2| that 
there is no mechanism to reduce the total energy of the system, therefore we look for the 
non-ground state condensation described by a macroscopic wave function that has nodes. 
We will classify the rotation curves for three thousand different visible matter distributions. 

II. QUANTUM FLUID OF DARK MATTER 

The system is described by a macroscopic wave function (order-parameter) ip. For sim- 
plicity, we consider spherically symmetric case. We also set the visible matter distribution 
to be spherically symmetric. This is not realistic, because the visible matter distribution of 
the disk galaxy is by definition not spherical. However, the purpose of this paper is to show 
the possibility of the decreasing curves and we believe that, for this purpose, it is enough to 
consider the spherically symmetric case only. We set the visible matter potential by hand 
and get the wave function of total system by numerical work. Also, we neglect inter-particle 
interactions apart from gravity. Since total distribution determines the potential of an in- 
dividual particle, we expect that the system be described by a solution of non-linear rather 
than linear wave equation. 

Newtonian potential V is given by 

V 2 V = 4%G(pdark + Pvisible) 

and the schrodinger equation is 

h 2 

ihd t ip = V 2 </> + V ■ ip(r) (1) 

2m 

By identifying p dark = GM m\ip\ 2 and defining p visi ue = GM mp v we get 

ihdttp = -— V 2 ^ + GmM I" dr\- f dr"4nr" 2 (\^\ 2 + p v ) ■ ip(r) (2) 
2m jo r l Jo 

Here, Mo is a mass parameter introduced for dimensional reason. We now set effective 
density of the visible matters as the Plummer potential for the bulge and exponential decay 
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for the disk. In eq. (1) it is important to remember that ip describes the state of the whole 
system rather than a state of an individual particle. The justification of the non-relativistic 
treatment can be found in the ref. 0. The basic idea is v ~ GmM ~ 10~ 2 , because in our 
case M ~ 10 12 M Q and m ~ 10~ 24 eV as we shall see. 

In this paper we consider two parameters that determine visible matter distribution. The 
one is the mass ratio of visible matter to dark matter and the other is the distance ratio of 
the distance in which 90% of visible matter exist to that in which 90% of dark matter exist. 

Now we come back to the Eq.(l). It is non-linear, hence we do not have the freedom to 
normalize the solution. However, as shown in ref. || this equation has a scaling symmetry 
which allows us to the resolve the ambiguity due to the choice of M . That is, the physical 
quantity M does not depend on the choice ofM . We are interested in a stationary solution 
xp(r, t) = e~ lEt / n ?p(r). We simplify the Eq.(l). After scaling by 

h 2 h 2 



r = r f , ip = r ip , p v = r p v , E = — e , r 



2m 2GM m 2 
We can write (1) in terms of the radial wave function u(f). 

u{r) + (e - f dr'lj f df"{u 2 {f") + p v (r")r" 2 )) ■ u{r) = (3) 



o r Jo 



where ip(r) = -^^r-- We solve this equation numerically. 



III. SCALING ANALYSIS 

In equation (2) there is a scaling symmetry. That is, the eq.(2) is invariant under the 
scaling 

u — > Am, f — > A _1 f, e — > A 2 e, p v — > A 4 ^. (4) 

Therefore, if u(f,e,p v ) is a solution, so is Xu(Xr, A 2 e, \ A p v ). Let M = M dark + M visih i e = 
j u 2 dr + J pr 2 . When we sacale f — ► A _1 f, M — ► AM. Since r = rgf is a physical 
quantity, it should be invariant under the scaling. Therefore tq — > A _1 ro- On the other hand, 
r ~ l/(M m 2 ), hence M m 2 —>■ XM Q m 2 . The scaling of the normalized velocity v can be 



found as follows. The virial theorem tell us v 2 = GM/r, so v = ^GM /r y M /r :— v • v 
Therefore under the scaling above, we should get t>o — > A _1 i>o to get physical velocity vqv. 
On the other hand, v ~ M m. Hence we get M m — > A _1 M m. Summarizing, under the 
scaling f — > A _1 f , 

M m 2 -> A _1 M m 2 
M m -> A _1 M m 

Therefore m is scaling invariant while Mo — > A _1 Mo. Consequently, M = Mq ■ M is scaling 
invariant as it should be. 

Why did we do above analysis? In solving Eq.(2), we must choose the two initial val- 
ues. The one is chosen u(0) = 0, because u(r) is a radial wave function. In case of linear 
Schodinger equation, there is a freedom of normalization so that one can freely choose first 
derivative of u at r = 0, or ^(0). Here, the equation is non-linear and we do not have the 
freedom of normalization. We are in trouble unless something happens. It is the scaling 
symmetry we proved that save our trouble. In our numerical work, whichever we choose Mo 
or A-Uo as the initial value, the physical quantity is not changed. We set u(h) = 0.05 * h, 
where h is a small segment of f. 

Now, we set the visible matter distribution as follows. 

P0 ' ((r/rA+l\* 5 f ° r illIler bul S e 

p v r = < u ' ; ; ' 

Pv(?d) ■ e^ rd ~ r ^ for the disk 

We study for various values of p an d fa- We classify the rotation curve (v(r)) with the mass 
and the distance ratio of the visible matter to the dark matter. 

IV. COMPARISON WITH OBSERVED DATA 

Now we present the main result of this paper. We looked for solutions with 5 or 6 nodes. 
This choice is motivated partly because, in this state there is a rotation curve very similar 
to NGC2998, and NGC801. Additional reason is if we take node number bigger than 7, 
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the total mass of the halo is larger than 1O 13 M . Hence there are some restrictions in our 
choice. Furthermore, in node number 4,3,2,1, the flat region is very narrow, and very hard 
to classify So, we choose node number 5, and 6. 

In figure- 1, we compare NGC2998 and 5.6M-19.7D where 5.6M means the ratio of (visible 
matter mass) / (total mass) is 5.6% and 19. 7D means (90% of visible matter dispersion 
distance ) / (90% of total matter dispersion distance) is 19.7%. There are four ripples in 
the NGC2998's data. But we suspect that outside the measured region there might still be 
density peaks and nodes that could lead to further ripples. We can see that the shapes of 
measured and computed rotation curves are very similar except near the core of the galaxy, 
where the error of measured data is relatively big. By comparing the maximal velocity and 
the size of a ripple, we can get v and r . By using m = h/(y/2r v ) and M = r v 2 /G, we 
get m = 2.90 x 10~ 2i eV and M = 7.117 x lO n M in this case. 

In figure-2, we compare NGC801 and 6.9M-12.7d. In this case m = 3.14 x 10~ 24 eV^ and 
M = 6.54 x 1O U M . 

In the figures-3 and 4, Y axis is visible matter's 'mass ratio' to total mass, and X axis 
is 90% of the visible matter 'dispersion distance' to the total distance. The figures 3, 4 are 
the classification of the shapes of rotation curves. The symbol / means that corresponding 
rotation curve is increasing type, similarly \ means decreasing, and - means flat. When we 
measure a rotation curve in a real galaxy, we can't see full shape since we can observe only 
to the extent that sufficient amount of the hydrogen gas exists. So up, down, flat (/, \ , - ) 
classification is not for the full shape, but for the region that corresponds to first four or five 
ripples. The figures show that in a galaxy with relatively big visible mass ratio and small 
extension ratio has decreasing rotation curves, while the one with small mass ratio and large 
extension ratio has increasing data. This result is entirely consistent with the observation 
that compact bright galaxies have decreasing rotation curves and small dwarf galaxies have 
increasing rotation curves ||. 

To give more intuition in reading the figures 3 and 4, we draw three rotation curves with 
the same mass ratio and different distance ratios. See figures 5, 6, 7, 8. We also give rotation 
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curves with the same distance ratio and different mass distance ratios. 

Notice that we can find some interesting fact on the visible matter distribution from 
figures 3, 4. Our analysis shows that the flat region is confined in a small window, mass 
ratio 2% to 10 % and distance ratio 10% to 20%. On the other hand, the overall flatness 
is a general phenomenon in the real galaxies. Consequently, in our frame work, we can say 
that most galaxies, which have flat rotation curves, should have the mass and distance ratio 
in this range. 

In figure 5, we fix the mass ratio 6.5% ±0.5%, and change the rate of dispersion distance. 
Figure 6 is node 6 case. As we change the mass fraction of the visible matter, we can find 
three different types of rotation curve. With the change of the mass fraction, we can obtain 
various types of rotation curves. If the visible matter distribution is highly confined, the 
shape of the graph is of a decreasing form. When the visible matter distribution is widely 
dispersed, the graph is of slightly increasing a form that is very similar to the case of dark 
matter only. When its distribution is neither very dispersed nor very condense, 

the graph is almost flat. 

In the case shown in figures 5,6, the mass fraction of visible matter is not so large. So, 
the shape of the graph is not so different to the case of dark matter only. However, when 
the mass fraction of visible matter is very large, the situation is quite different. See figures 
7, 8. In figure 7, we fix the mass ratio 25% ± 2%, and change the distance ratio in node 5. 
The same analysis for node 6 is shown in figure 8. Notice that in figures 7, 8 there is no flat 
curve in any case, and the shapes are quite different to the case of dark matter only. 

In figure 9, we fix the distance ratio 12% ± 1%, and change the mass ratio node 5 case. 
Figure 10, is for node 6 case. 

V. DISCUSSION AND CONCLUSION 

In this paper we investigated the rotation curves of galaxies assuming that the dark 
matter consists of ultra light boson appearing in late time phase transition. Galactic halos 
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made of this species are highly correlated bose liquid. We develop a Landau-Ginzberg type 
theory describing the collective behavior of the system. 

As the main results, we can obtain decreasing curves as well as flat and slightly increasing 
curves. In the previous work [0] we could see only slightly increasing curves, while observed 
data shows that a few galaxies have decreasing rotation curves. Here, we can see the flat 
curves and the decreasing curves as well by considering the visible matter, in addition to 
the dark matter. We show that the small amount of dark matter can change the shape of 
rotation curves near the core significantly. 

Our paper has one interesting prediction. The flat case is confined in a rather small 
window corresponding to the mass ratio 2% ~ 10% and in the distance ratio 10% ~ 20%. 
On the other hand of the observed galaxies has a 'flat rotation curve' So, we can predict 
that in most of galaxies have the mass ratio and distance ratio in that region. 

So we conclude that the ultra light boson appearing in the late time phase transition 
theory has good enough properties to be a dark matter candidate, at least from the rotation 
curve point of view. 

acknowledgement This work was supported by research fund of Hanyang University. 
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FIGURES 

FIG. 1. comparing the theory and observation for NGC2998 with the data that 5.6% of mass 
fraction and 19.7% of distance fraction 

FIG. 2. comparing the theory and observation for NGC801 with the data that 6.9%of mass 
ratio and 12.7% of distance ratio 

FIG. 3. classification of node 5 data; / means increasing, — means fiat, \ means decreasing 
rotation curve 

FIG. 4. classification of at node 6 data 

FIG. 5. (1) mass ratio = 6.4% and distance ratio = 28.0%. (2) mass ratio = 6.3% and distance 
ratio = 13.0%. (3) mass ratio = 6.2% and distance ratio = 5.9%. AT node 5 

FIG. 6. (1) mass ratio = 6.4% and distance ratio = 27.7%. (2) mass ratio = 6.3% and distance 
ratio = 13.6%. (3) mass ratio = 6.2% and distance ratio = 5.8%. At node 6 

FIG. 7. (1) mass ratio = 25.4% and distance ratio = 27.4%. (2) mass ratio = 26.8% and 
distance ratio = 13.8%. (3) mass ratio = 26.8% and distance ratio = 5.9%. At node 5 

FIG. 8. (1) mass ratio = 25.4% and distance ratio = 26.1%. (2) mass ratio = 25.8% and 
distance ratio = 14.3%. (3) mass ratio = 25.8% and distance ratio = 6.9%. At node 6 

FIG. 9. (1) mass ratio = 0%. (2) mass ratio = 8.4% and distance ratio = 12.1%. (3) mass 
ratio = 14.6% and distance ratio = 12.8%. At node 5 

FIG. 10. (1) mass ratio = 0%. (2) mass ratio = 8.4% and distance ratio = 12.2%. (3) mass 
ratio = 14.5% and distance ratio = 12.3%. At node 6 
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This figure "figl-6.png" is available in "png" format from: 
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